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A graphene quantum dot under intense ac field and static low magnetic field is investigated. 
From a tight-binding perspective, applying a Fourier-Floquet transformation and renormalization 
process, we observe that graphene -intrinsically anisotropic- reveals field polarization signatures in 
the quasi-density of states. For the ac field polarized along the armchair direction, the dressed 
electronic structure shows an emergent property: an ac field induced quantum ring. This is inferred 
by the orientation-dependent formation of a miniband of energy states periodically modulated with 
increasing magnetic field, exactly analogous to the behavior of a quantum ring spectrum. 
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An intense research effort has been launched after 
the stabilization of graphene^, a true two-dimensional 
system, due to the emergence of unconventional elec- 
tronic properties with broad potential applications^. A 
benchmark for a two dimensional system is the quantum 
Hall effect, which indeed was observed in graphene^, 
even at room temperatures^. The research effort on 
graphene also includes its response to intense ac electric 
fielda2»Si2ii£. The concurrence of an ac electric field in 
plane and a perpendicularly applied static magnetic field, 
acting simultaneously on a graphene system, will be the 
focus of our work here, motivated by another recently dis- 
covered emergent effect, namely the microwave induced 
zero resistance states in very high mobility, but still "con- 
ventional" , quantum Hall system s 11 ' 12 . One characteris- 
tic feature of these states is that they are immune to the 
in plane ac field polarizatio n 13 ' 14 , an aspect that con- 
tinues to be addressed theoretically in "conventional" 
2D system s 15 ' 16 . Furthermore, the microwave intensi- 
ties in these experiments could go beyond a perturbative 
regime, leading to the necessity of treating the electron- 
ac field interaction from a non-perturbative approach™. 

These elements settle the scenario and context of our 
work: we investigate the evolution of the graphene quasi- 
energy spectrum in the presence of both, an ac field paral- 
lel and a static magnetic field perpendicular to the carbon 
atoms plane. On the contrary to electrons at the bottom 
of the GaAs conduction band, graphene is intrinsically 
anisotropic and the response to ac fields should reveal 
fingerprints of this anisotropy. We address the prob- 
lem from a tight-binding approach applying a Floquct- 
Fourier transformation and renormalization procedure 
that permits us to find the density of states of the quasi 
energy spectrum^ 7 -. This approach leads to the limitation 
of handling the problem only for finite systems. On the 
other hand, emerging effects for dot like structures may 
be identified, and indeed, small quantum dots out of a 
graphene layer is rapidly becoming a realit y 18 ' 19 . 



The tight-binding approach for the graphene electronic 
structure coupled to an ac field is given by an Hamil- 
tonian divided in two parts, H = H a + Hi, where Hi 
includes the interaction of the lattice with the ac field, 
while the interaction of the static magnetic field with the 
lattice is included in H n , 
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where the index x is for the ac electric field compo- 
nent parallel to the zigzag direction, and y is for the ac 
field parallel to the armchair direction, as visualized in 
the bricklayer representation of a graphene finite lattice, 
Fig. [TJ Furthermore, (li,h) are the index of the (x, y) 
coordinates of the sites. The phase factor is defined as 
a = where $ e = h/e is the quantum magnetic flux 

and $ = a x a y B is the magnetic flux per a half unit cell of 
the graphene lattice, because a x = \/3a/2, in mean value 
a y = 3a/2, and a = 1.42 A, the graphene lattice param- 
eter. While the p z orbital energy will be taken constant, 
and we choose £i lt i 2 = with no diagonal disorder. The 
hopping parameter are V x and V y as the nearest neigh- 
bors interaction and are taken constant and equal to 2.97 
eV, along the continuous lines of the brickwall in Fig. [3^. 
The ac field is described by u> = 2irf and F, the frequency 
and the field intensity, respectively. 
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A Floquct-Fouricr transformation, considering the so 
called Floquet eigenstates |li,/ 2 , m ) = ^1^2) ® \ m ) with 
to indicating the photon number states, redefines the 
Hamiltonian H as an infinite matrix with the elements 
given bj*i£: 
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This matrix is truncated at dimensions given by 
L x L y (2M + 1). L x and Lj, are the lateral sizes of the 
graphene finite lattice in number of atomic sites, while 
M is the maximum photon index. Since the ac field 
couples the photon replicas, for example, the Floquet 
state indexed by m to states labeled by to — 1 or m + 1, 
multiple photon processes become relevant with increas- 
ing field intensity and M settles how many replicas are 
taken into account. This Floquet matrix is a tridiago- 
nal block matrix with L x x L y diagonal blocks, given by 
E M = (E — mhuj)\ + Ho representing a photon replica 
with the matrix elements given by the left hand side of 
Eq. [3J The coupling of system to the intense ac electric 
field is represented by the off-diagonal blocks J- with the 



elements given by T = e 



The prob- 



lem is numerically handled by means of a renormaliza- 
tion procedure^, based on the associated Green's func- 
tion. The final result is the dressed Green's function 
for one of the photon replicas, say to = 0, and a quasi- 
density of Floquet states, p{£ + in) can then be obtained 
as p(£ + in) = — A Im [ Tr Gmm ] , in the atomic sites 
basis. 



III. NUMERICAL RESULTS: LANDAU LEVEL 
REPLICAS AND FIELD INDUCED QUANTUM 
RINGS 

The construction of our system is based on the topo- 
logical equivalence between the brickwall-like and hon- 
eycomb lattices, as already pointed out by lye et ai21, 
and as illustrated in Fig. [T] for a rectangular dot of 7 x 6 
carbon-like sites. The actual system discussed in the fol- 
lowing is a rectangular dot of 21 x 20 sites^. The edges 
of this rectangular graphene dot are zigzag along the x 
direction with period a x and armchair along the y with 
period a y (see representation in Fig. []}. 

The starting point is the inspection of the density of 
states of the graphene quantum dot in the presence of 
a perpendicular magnetic field, but in absence of an ac 
field^S, shown in Fig. fS] In this figure both the electron 
(E > 0) and hole-like (E < 0) halfs of the spectrum are 



x axis 

FIG. 1: The bricklayer model for the graphene lattice. At the 
right(bottom) side an armchair(zigzag) edge is illustrated in 
the geometry of the real honeycomb lattice. 



shown. Along the horizontal axis at E=0 one can see the 
central Landau level (LL), n — 0. At higher(lower) ener- 
gies one clearly identifies the second, n = +(— )1, and the 
third, n = +(— )2, LLs with the typical y/nB dependence. 
Between the central and the second electron(holc)-likc 
LLs (n = +(— )1) we see the edge states going down(up) 
in energy with increasing magnetic field. At the up- 
per(lower) left corner, there is a wide region of strong 
mixing of the edge and bulk-like states, since there the 
magnetic length is comparable to the system size, (this 
is the so called weak field limit, as elegantly discussed 
previously for a square lattice model^). 

At higher magnetic fields, the right side of Fig. [21 the 
lattice effects start to develop giving rise to the self- 
similar Hofstadter spectrum for a honeycomb lattice. 
Our interest is focused on the energy versus magnetic 
field window in which the LLs are already well defined but 
the lattice effects are not yet relevant. Referring to Fig. [2J 
this window of interest spans from 0.02 < $/$ e < 0.04. 
Hence we are focusing on the energy and magnetic field 
window between the weak field limit and the threshold for 
the lattice effects. Indeed for the present small systems, 
bulk-like LLs start to be defined here at very large mag- 
netic fields. However, the present results can be scaled 
to lower magnetic fields in larger dots, since the proper 
length scale for defining a LL is simply a dot dimension 
larger than the magnetic length. 

The inclusion of an ac field with the electric field com- 
ponent in the graphene plane, parallel to either a zigzag 
or to an armchair direction is the main discussion stage in 
what follows. In Fig. [31 we show the electron-like part of 
the quasi density of states of the graphene dot electronic 
structure dressed by photons, defined by an energy of 
huj = 0.8 eV and a field intensity of 0.7 x 10 6 V/cm 
applied parallel to the zigzag direction. Initially these 
values have to be discussed in order to build up a clear 
perspective for the present work. This field intensity cor- 
responds to eaF f» 10 mcV, considering a in Eq.[3]as the 
lattice constant for graphene. These arc again very high 
energies and field intensities when compared to actual 
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FIG. 2: The graphene quantum dot spectrum as a function 
of a perpendicular magnetic flux. Color scale: the darker the 
line the higher the density of states. 



experimental conditions (the ac field are normally in the 
terahertz ranged for the effects under scrutiny). Never- 
theless, the high absolute values for the ac field intensity 
and photon energy also scale down for larger systems. In 
summary, the computational limitation to small systems 
lead the formation of LLs beyond the weak field limit to 
high magnetic fields and, in order to observe interesting 
photon dressing effects, the ac field parameters are also 
increased. Therefore the discussion of the results has to 
consider normalized units, E/hu and eaF/fko, leading to 
a view also valid for realistic and experimentally available 
parameter sets, as discussed below. 
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FIG. 3: The same spectrum of Fig. 2 now under an ac field 
applied parallel to the zigzag direction. The intensity of the 
field is 0.7 x 10 6 V/cm and the photon energy is 0.8 eV (/ = 
£;//i=193.44 THz). 



Reports on the microwave induced zero resistance 
states lead to estimates of the available field intensities 
in the 100 GHz range of eaF/Tioj f» 0.2. Such a ratio 
already points toward the necessity of non-perturbative 
approaches^. Hence, in Fig. [3] we limit ourselves to even 
lower intensities eaF/tko < 0.1, in order to carefully fol- 
low the modification induced by dressing with photons 
the bare electronic structure of graphene. In the spec- 
trum shown in Fig. [3] we clearly see the m = 1 replica of 
the n = LL (the horizontal line at E = Tiuj = 0.8 eV) 
while only faint signatures of the m = — 1 replica of the 
n = 1 LL can be perceived (leading to a shadowy region 
below n = 1 LL), as well as replicas of edge states. 




Magnetic Field (<t>/* e ) 

FIG. 4: Same as Fig[3] but with the ac field parallel to an 
armchair direction. 

Fig. Q] depicts the quasi density of states for an ac field 
with the same frequency and intensity, but now with the 
electric field polarized along the armchair direction. Al- 
beit the same field intensity, the replicas seem to be more 
intense, in particular the edge states ones, as can be par- 
ticularly seen in an energy stripe around 0.4 eV. This 
structure in the set of quasi edge states will be addressed 
in the following by means of further increasing the ampli- 
tude of the ac field. So far these results indicate a clear 
difference between the intensity of the photon replicas for 
fields applied either in the zig-zag and armchair direc- 
tions, but a clear description of the effect is still missing. 

In this context, the quasi density of states for the 
graphene dot under an ac field along the zigzag direction 
is shown for higher frequency and field intensity param- 
eters in Fig. [5] Tllo = 1.6 eV and eaF/hcu w 0.03. By 
comparing with Fig. 3 and having in mind the photon 
doubled frequency, two replicas of bulk LLs can be seen: 
now, besides the m = — 1 of LL n = 1, an analogous 
replica for LL n = 2 can also be seen. A very strong 
replica of the central LL can be seen at E = 1.6 eV, 
which now merges into the weak field limit states. 

Besides the anisotropy in the intensity of the replicas 
that could be experimentally identified, no clear coupling 



4 



between quasi states belonging to different replicas have 
been observed for this case of ac electric field parallel to 
the zigzag direction. Indeed, anticrossings between quasi 
states are expected in several systems leading to interest- 
ing effects like collapsing superlattice minibands due to 
intense ac field induced dynamic localization^ or, more 
recently and in the context of the present work, new gaps 
openings in graphenc systems^. Such couplings between 
different quasi states are revealed by increasing the inten- 
sity of an ac field polarized along the armchair direction, 
Fig. O Now the hole-like edge states associated to the 
m = 1 replica of the central LL (rising in energy with in- 
creasing magnetic field) anticross with the electron-like 
edge states associated to the bare (m = 0) LL (lower- 
ing in energy with increasing magnetic field). Therefore, 
these anticrossings occur at half the energy separation 
of both replicas, namely at hui/2 = 0.8 eV for the cho- 
sen parameters. This picture can be clearly identified 
by the gedankenimaging of superposing hole-like edge 
states shown in the lower half of the bare spectrum in 
Fig. 2 with the electron-like edge states of the very same 
bare spectrum. Such superposition, actually impossible 
in a bare system, becomes feasible with the building up 
of the LL replicas induced by the intense ac field. The 
appropriate anticrossing behavior leading to the forma- 
tion of a nicely defined band is connected to the proper 
symmetry of an ac field along the armchair direction. 

A careful inspection of the miniband build up by the 
anticrossings of hole-like and electron-like states reveals a 
periodic modulation with increasing magnetic field. The 
period of this modulation, A<f>/<f> e w 0.004, represents 
a flux quantum through an area corresponding approxi- 
mately 250 graphene unit cells, i.e., of the order of the 
dot size. This suggests that this miniband shows the be- 
havior of a quantum ring spectru m 20 i 26 near the edges of 
the dot. The simplest model for a quantum ring is a one 
dimensional tight-binding ring of sites, enclosing a mag- 
netic flux, which can actually be treated analyticall y 27 ' 28 . 
In experimental systems, like a nanostructured two di- 
mensional electronic gas in GaAs, quantum rings have fi- 
nite thickness but bona fide quantum ring spectrum can 
be identified with these very miniband modulations with 
magnetic field 2 ^. 

One should have in mind that such periodic modula- 
tion of the electronic structure is an outcome of the ring 
geometry, which is absent in our system. The actual po- 
tential felt by the electron is a photon dressed potential 
with a different effective symmetry^, here the case of 
a ring. From another point of view, the formation of 
a quantum ring like electronic structure without a ring 
structure but in the presence of an ac field is a unique 
feature of a graphene quantum dot. An ac field induced 
quasi quantum ring is only possible due to the hole-like 
edge states associated to the central LL m = 1 replica 
anticrossing with the m = electron-like edge states. It 
is worth to mention that it could be seen as an emergent 
property of the Dirac-like electronic structure, since such 
an effect cannot be observed for conventional two dimcn- 
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FIG. 5: A spectrum of a graphene quantum dot under an ac 
field intensity of 3.5 x 10 6 V/cm (eaF = 50 meV) applied in 
the zigzag direction of the lattice. Here the photon energy is 
1.6 eV (386.88 THz). 



sional electronic systems for which below the lowest LL 
there is only a magnetic barrier— and all LLs show a 
linear dispersion with magnetic field. 
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FIG. 6: The same as Fig. [5] with the ac field parallel to an 
armchair direction. We observe an ac field induced quantum 
ring spectrum of the coupled edge states. 



IV. DISCUSSION: THE QUEST FOR THE 
FIELD POLARIZATION EFFECT 



A remaining important question concerns the field po- 
larization dependence. At a first glance, like mentioned 
in the introduction, the consequences of the intrinsic 
anisotropy of grapheme could be explored particularly 
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in a low field quantum Hall regime, due to the insensi- 
tivity of microwave induced zero resistance states to the 
field polarization in GaAs based systems. 

Here qualitatively different pictures emerge with differ- 
ent ac field orientations. The quantum ring like spectrum 
only appears for the ac electric field parallel to the arm- 
chair direction. Nevertheless, the field induced LL repli- 
cas appear irrespective of the field polarization, albeit a 
quantitative difference in the intensity. Since the quan- 
tum ring like miniband is build up by means of proper an- 
ticrossings, rather involved symmetry properties should 
be involved. 

A further hint refers to the interaction between bulk 
LLs and their replicas. Looking to the region of the spec- 
trum in Fig. 6 where the m = 1 replica of the n = LL 
crosses with the bulk (m = 0) n = 1 LL, we see that this 
crossing occurs very near the weak field limit. Therefore, 
if we do not consider high photon energies, the investiga- 
tion of the coupling between these LLs is hindered by the 
presence of edge states. It is also due to the small size 
of the dot considered here (restricted by numerical costs) 
that we have to consider very high photon energies to 
follow the interactions between bulk LLs and their repli- 
cas. In spite of these features, the results suggest that 
the m = 1, n = LL replica seems to preserve its identity 
throughout the entire range down to the zero magnetic 
field limit, a very different situation from what is ex- 
pected for LLs with either a \JnB or linear dispersion 
with the magnetic fielcUt However, a closer observa- 
tion reveals th at this preservation of the identity of the 
LL replicas occurs only for the armchair case. For ac 
fields parallel to the zigzag direction, Fig. 5, the bulk LL 
merges with the whispering gallery of edge states in the 
weak field limit and is completely smeared out. An analo- 
gous smearing out happens to replicas of LLs with linear 



dispersion with the magnetic field, observed for square 
lattice models^. Therefore, electric fields parallel to the 
zigzag direction seem to have similar properties than elec- 
tric fields parallel to one of the sides of a square lattice. It 
should be mentioned that real material edges should be a 
mixture of zigzag and armchair edges and the robustness 
of the effect has to be addressed in further work on such 
ac field induced quantum rings. In actual geometrically 
defined quantum rings, the characteristic magnetic field 
periodic spectrum is robust considering edge disorder—. 



V. FINAL REMARKS 

The several results depicted here may be summarized 
in two main findings, (i) The unique electronic structure 
of graphene near the Dirac point, in the presence of a 
magnetic field, leads to anisotropic response to the ac 
field orientation, as revealed in the quasi density of states, 
(ii) The ac field induces a quasi-quantum ring (only for 
the ac field along the armchair direction) within the edge 
states spectrum, an emergent effect due to the flatness of 
the central LL and the special coupling among electronic- 
like and hole-like replicas of edge states. 
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